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Abstract. Data on multimode Cepheids from OGLE-III catalog of the LMC Cepheids are confronted with results from model 
calculations. Models whose radial mode periods are consistent with observation are not always in agreement with published 
evolutionary models. Nonradial mode interpretation is considered for the cases of unusual period ratios. The greatest challenge 
for stellar pulsation theory is explanation of double-mode pulsators with period ratios near 0.6. 
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INTRODUCTION 

The OGLE-III catalog of the LMC Cepheids [1] con- 
tains a large number of objects with more than one mode 
excited. There are 61 classical beat Cepheids pulsating 
in the fundamental and first overtone modes (F/IO) and 
203 objects pulsating in the first two overtones (10/20). 
The rare types of multimode radial pulsators include two 
10/30 objects, two triple-mode pulsators of the F/10/20 
type and three of the 10/20/30 type [2]. Of these rare 
Cepheids, only two 10/20/30 objects were known be- 
fore [3]. More common (29 objects) are first overtone 
pulsators with additional periods which are equal about 
0.6 of the dominant ones and cannot be explained in 
terms of radial overtones. We denote them as the lO/X 
pulsators. In Figure 1, we show period distribution of var- 
ious types of short period Cepheids. 

In a number of LMC Cepheids, additional periods 
close to the dominant modes are detected. These were 
found in a relatively small fraction (4%) of fundamen- 
tal mode Cepheids but in a significant fraction of ob- 
jects with exited overtones (e.g. 18% in lO and 36% in 
10/20). 

INFERENCE FROM DATA ON TWO 
RADIAL MODE PERIODS AND W/ 

Most of multimode pulsators are found among short- 
period Cepheids, which, as first noted in [4] and [5], are 
not well explained in terms of evolutionary models of 
core-helium burning stars (second and third crossings of 
the instability strip). Some of them may be in the first 
crossing phase, which is much faster, but there are also 
problems with this interpretation. 
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FIGURE 1. Period distribution for the LMC Cepheids with 
the first overtone excited. The 10/20 pulsators constitute 
13.8% of the sample but dominate in a certain period range. The 
F/IO pulsators constitute only 4.1%. The 10/X objects (1.9% 
of the sample) and the long-period F/IO pulsators occur in the 
same period range. 



Data on two periods of radial modes and on Wesenheit 
indexes, such as shown in Figures 2 and 3, yield strong 
constraints on stellar models. In principle, color data may 
also be used but they are less accurate. Instead, as a con- 
strain on effective temperatures, a requirement of mode 
instability was adopted in recent works exploiting data 
on multimode Cepheids ([6], [7], [8] and [9], hereafter 
DS). In DS, which was devoted almost exclusively to the 
10/20 LMC Cepheids, the Wi and period data (Figure 
3) were used to constrain parameters of these objects, in 
particular, their masses. 

In DS and here, model values of W] were obtained as- 
suming distance modulus to LMC 18.5 mag [10]. The er- 
ror of ±0.05 mag translates to ±7% uncertainty in stel- 
lar mass. Lines shown in Figures 2 and 3 were obtained 



0.76 


1 1 1 1 


ii | , 


1 , , , | , , , | i i 
3.0M S 3.5M 4.0M S 4.5M B - 






Aloj 


L = 0.0 - - 




□ 


Aloj 


L = 0.2 - 








L = 0.4 


0.74 




Aloj 


















-'■?i'.'! f » ... 


0.72 










, , 1 , 


, 1 , 


1 i i r I i i ' i 1 i i 





Alog L=0.0 
Alog L=0.2 
Alog L=0.4 
i I i 



-0.2 



FIGURE 2. Data (points) on the F/10 Cepheids. The Pe- 
tersen diagram (top) and PL relations (bottom) compared with 
values for selected models of indicated mass (segments). The 
segments correspond to the simultaneous instability of the F 
and lO modes for models of specified mass and luminosity. 
Luminosity increments AlogL refer to evolutionary models cal- 
culated with a ov = 0. 



for models calculated with the OP opacity data [11] 
for the new solar heavy element mixture [12], adopting 
Z = 0.006. We carried linear nonadiabatic calculations 
of radial modes for deep unfitted envelope models within 
the relevant temperature range using the code described 
in [13]. For luminosity, we first used the values obtained 
from evolutionary models in the first crossing phase, 
calculated assuming no overshooting. Then, we consid- 
ered models with luminosity increments AlogL = 0.2 
and 0.4 to reproduce the whole range of observed pa- 
rameters. Models with typical overshooting are brighter 
by AlogL ss 0.15. Those in the second and third cross- 
ings are brighter than models in the first crossing by 
AlogL w 0.25 (see [14]). 

In the lower panels of Figures 2 and 3, we may see that 
models of constant mass form nearly a single line in the 
logf — Wj plane for specified mode and metallicity. In 
fact, the dependence on metallicity is quite weak. This 
means that, once we have credible Wj and the LMC 
distance, we have quite accurate assessment of a Cepheid 
mass. 

The upper panel of Figure 2 shows the classical Pe- 
tersen diagram. We may see that the measured period 
ratios are well reproduced with our models. Compar- 
ing plots in the two panels, we note certain problem for 
logPp < 0.3. Models with masses near 3M Q are in agree- 
ment with the observational PL relation if some luminos- 



FIGURE 3. The same as Figure 2 but for the 10/20 
Cepheids. 



ity excess (presumably due to overshooting) is allowed. 
Thus, the objects may be understood as first crossers. The 
calculated period ratio is somewhat low. The agreement 
would improve by choosing lower Z value (see [8]). For 
the F/IO Cepheids with logi^i > 0.3, both the period ra- 
tios and the Wesenheit indexes, point to masses near4M 
and the luminosity excess AlogL w 0.4. These are pre- 
sumably helium burning objects. 

Figure 3 is from DS, where it was concluded that 
great majority of the LMC 10/20 Cepheids have masses 
M = 3.O±O.5M . The objects with logPio < -0.2 are 
well explained with post-main sequence stellar models 
crossing the instability strip for the first time, calculated 
assuming no or moderate overshooting. For those with 
longer periods, which constitute the majority of the sam- 
ple, a significant luminosity excess is needed. Interpre- 
tation in terms of first crossing models requires an over- 
large overshooting, perhaps connected with fast rotation. 
The difficulty of this explanation is the short crossing 
time of the instability strip while the number of the ob- 
jects showing the luminosity excess is relatively high. If 
these are helium burning objects then the difficulty is the 
low inferred mass. Standard evolutionary tracks for stars 
with M < 3.5 and acceptable metallicity do not enter the 
instability strip in this evolutionary phase. 

All five triple-mode Cepheids are most likely first 
crossing objects. Detailed modeling of two 10/20/30 
objects implied masses near 3.3M [6]. The third one has 
very similar periods and Wesenheit index, hence, must 
have a similar mass. Models of similar mass and calcu- 
lated with moderate overshooting well reproduce data on 



the two 10/30 Cepheids. The two F/10/20 objects oc- 
cur at very different periods. The one at logPjo = —0.4 
has mass M < 3.OM0 and that at logPio = 0.03 has mass 
M«3.5M . 



NONRADIAL MODES AND 
TROUBLESOME PERIOD RATIOS 

In the radiative interiors of evolved stars, the Brunt- 
Vaisala frequency, N, is by orders of magnitude larger 
than frequencies, CO, of radial modes. This implies that 
all nonradial modes propagate as high-order gravity 
waves with radial wave number, k r , satisfying rk r — 
y/£(£+l)~ 3> 1. The displacement amplitude is | 
[C+e^ + C-e-'^e- 10 " with ^ = Jdrk r . 

Large k r leads to large radiative losses. Global modes 
exist as long as |C_| ~ \C+\, and they are unstable if the 
driving effects in outer layers are strong enough. How- 
ever, the wave braking occurring at k r \^ \ ss 1 may pre- 
vent growth of an unstable mode amplitude to observ- 
able level. Once it occurs, the reflected wave must be ig- 
nored (|C_| « 0) and this implies wave energy losses. 
Only with a resonant energy transfer from a nearby ra- 
dial mode, detectable amplitude of nonradial modes may 
be maintained. This mechanism, suggested as an expla- 
nation for peaks close to radial mode frequencies in RR 
Lyrae stars [15], may apply also to Cepheids. 

Strong radiative damping in the interior leads to 
C_| <C \C+\. This again justifies ignoring the reflected 
wave. At sufficiently high-£, there are well-trapped 
modes, which, in spite of the wave losses, remain unsta- 
ble [16], [17]. Due to cancelation of opposite sign contri- 
butions, the net light changes resulting from such mode 
excitation are expected low. 

The Px/P\o period ratios are concentrated around 0.6 
and 0.62. There are no systematic differences in the 
positions in the PL diagram between the two groups of 
the IO/X Cepheids. Both are well inside the IO band. 
No radial mode may account for such period ratios. The 
closest but still far away is P30/P10, as we may see in 
Figure 4. 

Interpretation in terms of nonradial modes is also dif- 
ficult. In models constrained by observational data, the 
only modes within the Px range which do not suffer 
from strong damping in the deep interior are F-modes (no 
nodes outside g-wave propagation zone) of i between 40 
and 50. For modes of degrees I > 6, the light changes 
arise mainly from geometrical distortion of the surface. 
Then, for even degree modes, the net light amplitude de- 
creases with I only as £~° ,5 . Still the amplitude reduction 
is large and would require intrinsic pulsation amplitudes 
of these modes significantly larger than those of the IO 
modes. Moreover, we have no explanation for the prefer- 
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FIGURE 4. Period ratio (top) and Wesenheit indexes (bot- 
tom) for IO/X Cepheids (symbols) compared with calculated 
values for the 10/30 Cepheid models (segments). Unfitted en- 
velope models were calculated for temperatures corresponding 
to the 10 instability range (30 is stable in these models). 



ence to 0.6 and 0.62 period ratios. 
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